Harmonic Production in Ferromagnetic Materials at Low 
Frequencies and Low Flux Densities 

By EUGENE PETERSON 

Synopsis: When a multi-channel communication circuit includes a non- 
linear element sucn as a ferromagnetic core coil, distortion of the wave form 
impressed upon the circuit is produced. In terms of the single frequency 
components, this distortion is manifested in the appearance of new com- 
ponents. This distortion may give rise to the reduction of quality in any 
channel, and it may also introduce crosstalk and interference, which consists 
of new frequencies not present in the impressed wave of any channel under 
consideration, produced by independent channels. In view of the recent 
increased use of multi-channel systems, it has become necessary to in- 
vestigate the effects of this type of distortion, to determine the dependence 
of this distortion upon the properties of the magnetic materials constituting 
the cores of inductance coils and transformers, as well as upon the circuit 
impedances, and to determine those constants of core materials which are 
significant in the distorting process. 

The behavior of magnetic materials to complex waves of magnetizing 
force is ordinarily a highly involved process, so that a direct correlation 
between distortion and some of the easily measured constants of materials 
is a matter of some difficulty. It has been established experimentally, as a 
confirmation of theoretical speculations, that the third harmonic e.m.f. 
generated by a sinusoidal wave of magnetizing force may serve as an index 
of the distortion with a complex wave of magnetizing foice. This relation 
is valid for low flux densities and for frequencies at which the screening effect 
of eddy currents is not important. The paper is therefore devoted to an 
investigation of the third harmonic production in its dependence upon the 
properties of hysteresis loops. These loop constants in turn are shown to be 
deducible from AC bridge measurements on a coil of known dimensions hav- 
ing a core of the magnetic material under investigation. The loop con- 
stants for a few materials are included in the text. An analogy exists be- 
tween the treatments of hysteresis loop and of three-element vacuum tube 
characteristics which enables us to compare simplifying relations introduced 
by Rayleigh and by H. J. van der Bijl in the two cases. 

The theoretical deductions are found to be in general agreement with 
experiment, and are applied to a number of cases of practical interest. 
These include the effects of air gaps and dilution, and the choice of core 
material in third harmonic production by inductance coils and transformers. 
Finally, the amount of third harmonic current flowing out of long lines is 
deduced with both lumped and continuous loading. 

Part 1. Hysteresis Loops and their Mathematical 
Representation 

NEW and improved systems of multi-channel communication which 
have come into use during the past few years have imposed 
rigorous requirements on the circuit elements constituting the com- 
municating link, and have made it necessary to investigate the degree 
of distortion which arises from the use of ferromagnetic apparatus. 
The distortion introduced may have two general effects: distortion 
of the signal in any one channel, which is usually the minor effect, 
and production of crosstalk and interference between the various 

762 



HARMONIC PRODUCTION IN MAGNETIC MATERIALS 763 

channels to which the non-linear element is common. Such elements 
are, in general, ferromagnetic core coils or transformers. The distor- 
tion introduced by the non-linear relation between flux density and 
magnetizing force is therefore of fundamental importance in the design 
of iron core coils and transformers which carry simultaneously a 
number of communication channels. 

The use of iron core coils and transformers in communication work 
is confined to comparatively low flux densities in contrast to the ordi- 
nary practice in power work, where operation usually occurs above the 
knee of the normal magnetization curve at a value of the order of a 
thousand times greater than that used in communication work. There 
are two main reasons for this restriction: losses are reduced, and the 
relation between flux density and magnetizing force approaches line- 
arity so that distortion is minimized. 

Under actual operating conditions in which the more important 
crosstalk effects arise, we have a complex wave of magnetizing force 
acting on the magnetic core. Non-linearity in the magnetic circuit 
gives rise to new frequencies which normally 1 are related to those 
impressed upon the circuit, being sums and differences of integral 
multiples of the originally impressed frequencies. These modulation 
products are all of odd order - when the core is unpolarized. On 
purely theoretical grounds we would expect to find relations between 
the amplitudes of the different frequencies resulting from any one 
order of modulation. To take the third order modulation products of 
two impressed frequencies (f u / 2 ) as an example, we would expect the 
amplitudes of the harmonics 3/\ and 3f-> to be related to the amplitudes 
of the other third order products: 2/i ±/o, 2/ 2 ±/i- This has been 
confirmed by direct experimental test, so that we are enabled to use 
the generated third harmonic voltage as an index of the generated 
voltages corresponding to the other third order products. Accordingly, 
we shall deal in the following with the third harmonic produced by a 
sinusoidal wave of magnetizing force, and so avoid a more involved 
analysis. 

The fundamental relation in the operation of ferromagnetic appa- 
ratus is of course the relation between the flux density B and the 
magnetizing force //. In contrast to the usual behavior of circuit 
elements, the relation between the two fundamental quantities — the 
independent and the dependent variables, II and B in this case — is a 
function not only of the value of the independent variable, but is also 

1 This is true in the low flux density region. At high densities and with highly 
reactive circuits as in magnetic modulators, other frequencies are sometimes found 
which correspond to natural oscillations of the coil and circuit. 

- Bell System Technical Journal, Jan. 1928, pp. 110, 111. 
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a function of its previous history as manifested in the phenomenon 
of hysteresis. This complicates matters to an extent far greater than 
is the case with other circuit elements, and some analysis has been 
carried out in which the hysteresis loop has been replaced by the 
normal magnetization curve, or by some such single valued relation 
between the variables. For some purposes this convenient simplifica- 
tion — it cannot be called a close approximation for our present pur- 
poses — is satisfactory, while for others it does not begin to tell the story. 
Inasmuch as our aim here is to deal with the actual phenomena 
involved rather than to arrive at some arbitrary procedure for repre- 
senting the facts, we shall in the following base considerations upon 
the hysteresis loop. 

Fig. 1 will serve to illustrate the effect of previous history upon 
flux density. The main loop there, which extends between the two 




Fig. 1 

limits of magnetizing force — H and H, is obtained when the magne- 
tizing force is varied cyclically between those two values in such a 
way that the magnetizing force has but one maximum and one mini- 
mum per cycle. In that case the B-II loop is traversed in the direction 
shown by the arrow. It is independent of frequency when the eddy- 
current losses in the iron are small, as we shall suppose them to be, 
and it is independent of the wave form of the magnetizing force so 
long as that wave form satisfies the condition we have laid down above. 
A sinusoidal magnetizing force, for example, satisfies that condition. 
When the magnetizing force contains components of such magnitude 
and phase that the wave form has multiple maxima or minima, the 
simple B-H loop no longer suffices to represent that relation, but 
auxiliary loops shown in the figure are involved. 
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Thus suppose a subsidiary maximum to exist at hi', as the magnetiz- 
ing force decreases, the flux density no longer follows the main loop 
but branches off on a subsidiary loop as indicated by the arrows. 
When the subsidiary minimum at ho is reached a new branch is started 
which completes the subsidiary loop, and which brings the magnetizing 
force back to the main loop at the point from which it originally 
diverged, and the main loop is thereafter followed until another 
maximum (or minimum as the case may be) of the magnetizing force 
wave is reached. For simplicity in the following we are going to deal 
solely with a sinusoidal wave of magnetizing force so that subsidiary 
loops are never called into play. With this understood, the relation 
between B and H is described by the simple loops of Fig. 2 over a 




Fig. 2 



certain range of magnetizing force, each loop being defined by a 
particular value of maximum magnetizing force. Each loop may be 
considered as constituted by two branches which join at the maximum 
field of the loop. 

It is clear, therefore, that with a periodic magnetizing force having 
but one maximum and one minimum per cycle, the flux density depends 
upon three properties of the magnetizing force — the maximum value, 
the instantaneous value, and the sign of dh/dt, being located on the 
lower branch when dh/dt is positive, and on the upper when it is 
negative. Now it is of course evident that when a definite loop form 
is available, a numerical solution by graphical or step-by-step methods 
may be had. It is further evident that, in the case of a definite im- 
pressed magnetizing force, the B-H loop may be broken up into its 
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harmonic components, as was done by S. P. Thompson. 3 Solutions 
in this form possess all the advantages and disadvantages of numerical 
ones in which any change of conditions leads to a new problem; the 
solution desired is a general one which will describe the phenomena in 
terms of coefficients characteristic of the magnetic material, and this 
type of solution is the subject of analysis in the following pages. 

Perhaps the least difficult method of arriving at an analytical 
solution without making any assumptions as to the form of the loops 
is the following. A power series for each branch of any loop is 
formulated, and the two resultant equations are combined in a trigo- 
nometric series. In this way the solutions, each one valid over but half 
the cycle, are combined to represent the relation of B to II over the 
entire cycle. 

The flux density on each branch of the loop may be defined in 
terms of two values of magnetizing force, one the maximum value of 
magnetizing force on the particular branch with which we happen to 
be concerned and the other the instantaneous value of the magnetizing 
force on that branch. The equation for either branch may then be 
expressed as a double power series in these two variables, the instan- 
taneous magnetizing force (h) which will be expressed in gilberts/cm, 
and the maximum value of the magnetizing force (77), expressed in 
the same units; 

00 00 

B(h,H) = £ £> mB MT», (1) 

7/1=0 n=0 

where 

1 dB(h, H) 



a mn = 



mini dh m dH n 



(2) 



The parallelism of this representation with that for the plate current- 
grid potential curves of a three electrode thermionic tube is evident, — 
in both cases double power series are involved. 4 The coefficients a m „ 
are derivatives which are evaluated at the point // = 0, II = 0, and 
it will be understood in (2) that these particular values are inserted 
after the derivatives have been taken, in quite the usual manner. 

There is a further interesting parallel here to the vacuum tube case 
regarding simplification of the general relation. In the case of the 
vacuum tube a simplification of the double power series due to van 
der Bijl has been employed which represents the family of tube char- 
acteristics by an equation in a single variable. This simplification 
consists in assuming the amplification factor to be constant, and the 
important point for us here is that it is equivalent to the assumption 

3 "On Hysteresis Loops and Lissajous' Figures," Phil. Mag., 1910. 
1 Peterson and Evans, "Modulation in Vacuum Tubes used as Amplifiers," Bell 
System Technical Journal, July, 1927. 
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that the different branches of a family have the same form, so that by 
suitable change of plate or of grid potential the plate current-grid 
voltage curves may be superposed. A relation of precisely the same 
form in the case of magnetic hysteresis loops was stated by Lord Ray- 
leigh 5 upon examination of data obtained by a magnetometric study 
of the behavior of a single low permeability specimen. Examination 
of his data enabled Rayleigh to conclude that the branches corre- 
sponding to different loops of a family could be superposed when re- 
ferred to a common loop tip, the branches all having the same parabolic 
form within the limits of accuracy of the measurements, over the 
range of magnetizing forces involved. 

It is of course evident that even if this relation held at low fields in 
all materials it would break down at sufficiently high fields. Further, 
there is no a priori reason to expect this relation to hold for magnetic 
materials other than the one Lord Rayleigh investigated unless we 
restrict consideration to a very small range of magnetizing force. 
With these ideas in mind it seems the safer procedure to assume no 
such simple relation between the different loops of a family, however 
convenient it might be, and to treat the problem in more general terms ; 
if any such simplifying relations exist they will be made apparent 
after application to definite materials. We may anticipate matters a 
bit to state at this point that certain materials seem to obey the 
relation while certain others seem to violate it within the range of 
forces involved in communication work, and further light is shed on 
the significant processes involved, in harmonic production by treating 
the problem in this way. 6 

General Equations for Hysteresis Branches. The equations of both 
the upper branch family and the lower branch family have the form of 
equation (1) — we may designate the upper branches by Bi and the 
lower branches by B 2 — but the coefficients of the two series differ in 
general. 

In order to put the equations in shape so that they may be of 
practical utility it is now necessary to determine the coefficients of 
the expressions for B x and B-> so that they apply to a definite loop 
family, and this is accomplished by reference to some of the more 
general properties of the loops and of the normal magnetization curve. 
These properties are as follows: 

6 "Notes on Electricity and Magnetism, III," Phil. Mag., 1887, V. 23. 

6 Unpublished work based on assumptions which include Rayleigh's relation was 
independently carried out by W. P. Mason of these Laboratories in 1922 and 1923. 
An account of some applications of Rayleigh's relation is to be found in an interesting 
paper by Jordan published in the Elektrische Nachrichten Technik, B. 1, H. 1, July 
1924. 

49 



768 BELL SYSTEM TECHNICAL JOURNAL 

1. When both h and H are zero the flux density on either branch 

is zero. 

2. The flux density on one branch with II and h given is equal and 

opposite in sign to the flux density on the other branch corre- 
sponding to the negative of h and to the same maximum 
force //. 

3. The two branches corresponding to a definite H meet at the 

normal magnetization curve. 

The application of these properties to the power series enables us 
to deduce relations between the coefficients as demonstrated in 
Appendix 1 : 

floi — floo = 0, 

002 = — 020, (7) 

003 = — 021- 

We shall find it sufficient to include the third degree terms for our 
work, so that we need not investigate relations between coefficients 
of higher degree. The loop equations are simplified by utilizing (7) 
and we can make a number of interesting deductions. Thus the equa- 
tion for the normal magnetization curve is given by 

B(H, H) = a 1Q H + a u H 2 + (a 12 + a 30 )H 3 . (6a) 

In this equation flio will be recognized as the initial permeability 
usually expressed as fi since upon division by II we have for the 
permeability 

M = 0io + a u H + (012 + a 30 )H 2 

According to this equation the change of permeability with magnetizing 
force is linear at sufficiently small fields. The above equations are, 
in all rigor, infinite series, but for our purposes it will be found sufficient 
to consider only coefficients of the third and lower orders, — in some 
cases the second order will suffice. 

An expression for the remanence curve of the loop family may be 
obtained by setting h equal to zero in the equation for the upper 
branch. In that case we find from (4a) 

B(0, H) = a 02 H°- + aosIP + • • • , (8) 

hence for sufficiently small magnetizing forces the remanence increases 
as the square of the magnetizing force. 

The hysteresis loss per cycle per unit volume may also be obtained 
directly from the branch equations (4a) and (5a). The loss in ergs, w, 
is equal to the area of the hysteresis loop divided by 4x. If now we 
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consider the loop area to be built up of strips of infinitesimal width 
based on the /j-axis, the height of any one of the strips is given as 

B = B x (h, H) - Bo(h, H) 
and we have 

1 r H 



1 r" 
= f- | Bdh (9) 



so that, by Appendix 1, (9) may now be expressed as 

w = £- (aoiH 3 + a Q JP +•••)• (10) 

It is clear, therefore, that at sufficiently low fields the hysteresis loss 
varies as the cube of the magnetizing force, and diverges when the 
field is made sufficiently large — it seems to be in general agreement 
with experimental results on ballistic loops and, as will be pointed out 
later, is verified by impedance change data under alternating excita- 
tion. In view of the remanence curve equation, (10) may be rewritten 
as 

w=^-IIB(0,H). (11) 

Various approximations have been made in the past to hysteresis loop 
forms in order to obtain convenient expressions for the hysteresis loss. 
Thus if we consider the loop as an ellipse the loss becomes 

\hb{o,h), 

while if we consider the loop a parallelogram the loss is 

- HB(0, H). 

Both these expressions give too large a result since the coefficient 2j3nr 
of the exact equation (11) is 0.212. 

Branch Equations for Materials Obeying Rayleigtis Relation. Ray- 
leigh's observation enables us to establish relations between the 
different coefficients involved in our development above when the 
loops of a family are similar in form. For the sake of completeness a 
derivation of these relations is given in Appendix 2 ; their validity 
may then be judged by test in specific cases. In the derivation we 
assume a power series expansion for one branch of the largest loop 
of a family referred to the tip, and assume that the smaller loops 
are of the same form. Then by referring the equations to the origin 
instead of to the loop tip, we arrive at the hysteresis branch equations. 
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By comparing the coefficients of this equation with those of the general 
equation (1), we can deduce relations between the coefficients. These 
are, from Equation 16 of Appendix 2, 

aoi = 0, 

#20 + #02 = 0, an = 2a2o, /4^\ 

#21 + #03 = 0, ai2 = #2i = 3aao, 

#40 + #22 + #04 = 0, an = #13 = 4a 4 o = 6a 2 2- 

The left column is in agreement with Equation (7), while the right 
column furnishes new relations between the coefficients. Thus the 
coefficients based on loop similarity are not inconsistent with those 
deduced under no assumptions, but the former involve additional 
relationships between coefficients which need not be satisfied in the 
general case. 

Families of hysteresis loops obtained by the usual ballistic method 
have been examined for these relationships in the case of several 
materials with the following results for coefficients up to and including 

the second degree : 

Hysteresis Loop Coefficients 





Silicon Steel 


'B' Dust 7 


•C* Dust 7 


Ol0 
Oil 

an/2 

C20 


270 

2600 

1300 

967 

- 951 


35.3 

1.53 

0.76 

0.65 

-0.71 


26.2 

0.59 

0.29 6 

0.29 

- 0.28 5 



The values tabulated were obtained by 'smoothing' data obtained 
from ballistic loops, which leaves a great deal to be desired from the 
standpoint of precision. The difference between a 02 and — fl 2 o is an 
index of this lack of precision. The average deviation of an/2 from 
the mean of a z and 020 (a relation deduced on the basis of similarity) 
is small for C dust and large for silicon steel while the third order 
coefficients show much poorer results. The lack of precision in the 
original data however, is such as to leave unsettled the question of 
loop similarity for the materials tested at the fields to which the 
coefficients apply ; the experimental error is too great. 

The use of ballistic methods for determining coefficients becomes 
even less satisfactory for materials with low hysteresis losses and 
cannot be used for analyses which have any aspirations to precision. 
It will be shown in the next section, however, that the significant 
coefficients enter into the impedance of a coil which employs the 
material under examination as a core, so that the desired coefficients 

7 Speed and Elmen, "Magnetic Properties of Compressed Powdered Iron," 
A. I.E. E., V. 40. 
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may be obtained with the aid of AC bridge measurements under 
appropriate conditions for which eddy currents and winding capacities 
are unimportant. These measurements are of a higher order of 
precision than those obtained by the ballistic method, and will be 
treated in detail in the next section. 

Part 2. Alternating Magnetization 
Sinusoidal Magnetizing Force. The same magnetic characteristics 
which were the subject of the investigation of the last section are 
involved in alternating magnetization when the applied field has but 
one maximum and but one minimum per cycle, and when the eddy 
losses are not great enough to introduce screening effects. As one of 
the results of that investigation, we have arrived at equations for both 
the upper and the lower branch families, so that to obtain an expression 
for the flux density valid over the entire cycle, it is now necessary to 
combine the two branches by a Fourier's series in the usual manner. 
Before doing this, however, it is necessary to express the branch 
equations for the flux density in terms of time, rather than as series of 
powers of the independent variable, the magnetizing force. To do 
this we substitute the equation for the magnetizing force 8 

h = 11 cos pt (18) 

in the branch equations (4a) and (5a) of Part 1. The result of the 
substitution is to give the upper and lower branch flux equations in 
terms of powers of cos pt, which may be expressed in terms of multiple 
angles. These two equations are then combined in a Fourier series 
which is valid over the entire cycle: 

B = ?? + E (°k cos kpt + a k sin kpt), 
z 1=1 

in which the coefficients are given by the usual expressions, equation 

22a of Appendix 3. The coefficients b , b 2 , • • • b 2k , a , a 2 , • • • a 2k are 

found to be identically zero on account of the symmetry of the loop 

family about the origin, while the fundamental and third harmonic 

coefficients are found to be as follows: 

ai = ^- (a 02 IP + aozH 3 ), 

b l = aioH + a u ff 2 + (a 12 + $a 30 )IP, (25) 

a 3 = - Tc~ (a<>2# 2 + aoalP), 
15ir 

b 3 = a 30 H 3 /4. 

8 The purely sinusoidal magnetizing force may be obtained, despite the varying 
reaction of the iron core coil, by connecting a generator through a low pass or band 
pass filter having a high impedance outside the pass band, or by connecting a pure 
sine wave generator to the coil through a high impedance. 
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Details of the derivation are given in Appendix 3. Inasmuch as we 
assume the applied field to vary as cos pt, the b's are in phase with the 
applied field and the a's are in quadrature. The two a's, it is observed, 
are connected by a constant of proportionality (a x = — 5a 3 ) and 
depend upon the remanence; or upon what is the same thing, the 
hysteresis loss divided by the magnetizing force. 

If we expanded the loop equations to higher powers we should find 
that <2i and a 3 cease to be linearly proportional. The coefficient &i 
is observed to have its first three terms identical with those of the 
normal magnetization curve, but the fourth term differs by precisely 
the amount b 3 . 

The voltage existing across a coil enclosing a core characterized 
by the above coefficients is 

E = nAlO- 8 ^- , (26) 

where n is the number of turns enclosing the core, A is the core area 
in cm 2 , B is the total flux density and E is the total generated potential. 
Carrying out this operation we have 

E = nAlO~ s (pai cos pt + 3pa 3 cos 3pt 

— pb\ sin pt — 3pb 3 sin 3pt) (26a) 

in which the voltage components in phase with the current depend 
on the hysteresis coefficients, and the quadrature components depend 
only on the coefficients of the normal magnetization curve. The 
dependence of each of these components upon the applied field is clear 
from Equation (25). To take the two third harmonic components it 
will be observed that a 3 starts to vary with the square, while b 3 starts 
to vary with the cube, of the applied field. It follows therefore that 
at sufficiently small amplitudes the third harmonic is produced by the 
a 3 term and not by the b 3 term, which means that under the conditions 
noted, harmonic production is due to hysteresis and not primarily to 
permeability change. 

From the two fundamental components of voltage across the coil 
an expression for the inductance and resistance offered to the flow of 
alternating current may be deduced. The inductance, it is easy to 
see, is obtained directly from the magnetization curve, — the d.c. 
permeability of the normal magnetization curve therefore coincides 
with the a.c. permeability at small fields. The resistance may be 
obtained from the expression derived above for hysteresis loss per 
cycle per cc. If we multiply that value by the volume of iron in the 
coil, by the frequency, and by 10~ 7 to convert ergs to watts we have 
the loss per second, and this may be equated to the square of the 
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effective fundamental current multiplied by the hysteresis resistance, or 

wfAwdlO- 7 = PR,J2. 

Here d is the diameter in cm, i" is the peak value of the fundamental 
current, / is the frequency, and R h is the hysteresis resistance. This 
may be solved for the hysteresis resistance in terms of the r.m.s. 
value of fundamental current /, which is at low fields, 

R h = 1.2-10r*n'Af7 am /d*. (27) 

The hysteresis resistance at small fields thus varies linearly with the 
applied current. 

It is easy to derive the relation between the hysteresis resistance 
and the generated third harmonic voltage at low fields since they 
involve the same constants a 02 and a Q3 ; from (26a) and (25) we have 
for the r.m.s. third harmonic voltage 

E 3 = 3pnA10- s a 3 /^ = 0.72 ■ 10-*fn 3 Al*a 02 !d 2 . 

Comparison with (27) shows that 

E 3 = 0.6RJ. (28) 

This simple relation is valuable in obtaining an idea of the harmonic 
production in a specific coil through resistance measurements, and 
more than that, it enables us to determine the coefficients significant 
in the distorting process for the material under test, so that the har- 
monic production in a coil of any dimensions enclosing the same core 
material may be calculated. The degree of precision ordinarily attain- 
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able is brought out for two materials in Fig. 3, and the agreement is 
observed to be within the experimental error in those two representa- 
tive cases. 
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These relations provide us with an expression for the ratio an/a 02 
which may be obtained from impedance measurements. It will be 
recalled that Rayleigh's relation would give this ratio the value two, 
and a.c. bridge measurements may be used to check this relation. 
The change of resistance with current is given by (27), and the change 
of reactance with current may be calculated from (25) and (26), since 
we have 

AX = a n H 2 pt7A10-*/l-y!2 

= lA2-10- s n 3 AfIa n /d 2 . 
Hence if we denote Rh by AR, we may write 

fl02 Aic 

The results obtained on some dust cores of different composition and 
two solid materials may be tabulated as follows: 

Material Number of Specimens On/ao2 

Grade ' B ' Iron Dust 10 2.65 

Grade ' C ' Iron Dust 1 2.63 

Permalloy Dust 9 'A' 3 2.10 

'B' 3 1.81 

Iron Dust No. 1 1 2.80 

Iron Dust No. 2 2 2.55 

Perminvar 10 1 3.0 

Alloy 'D* 2 1.84 

The method of analysis and the results obtained above may be 
summarized as follows. Starting with the general development of 
the hysteresis branches in a double power series and making use of 
only the most fundamental experimental observations, equations have 
been derived for the normal magnetization curve, the remanence 
characteristic, and the hysteresis loss characteristic in terms of the 
original hysteresis branch coefficients by combining the equations 
for the two branches in a trigonometric series. The series for the 
normal magnetization curve is found to start with the first power of 
the magnetizing force, the remanence starts with the second power, and 
the hysteresis loss starts with the third power, the curves varying in 
accordance with their respective first terms when 7/ is small. These 
results seem to be in general agreement with experience. For large 
values of II the shapes of the different curves depend of course upon 
the size and sign of the coefficients involved, about which nothing 
can be said until the coefficients have been evaluated from experimental 
data. 

9 Shackelton and Barber, "Compressed Powdered Permalloy," A. I. E. E. Con- 
vention, February 1928. 

10 Elmen, "Magnetic Properties of Perminvar," /. F. I., September 1928. 
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In the case of a sinusoidal magnetizing force the fundamental and 
third harmonic flux densities were derived with the aid of the trigono- 
metric series development. These were shown to depend in a simple 
way at low fields upon the normal magnetization and hysteresis loss 
characteristics. The in-phase fundamental voltage component de- 
pends directly upon the hysteresis loss per unit current, while the 
fundamental component in quadrature varies with the magnetization 
curve at low forces. The in-phase third harmonic varies with the 
in-phase fundamental at low forces, and the third harmonic in quadra- 
ture comes in only with larger forces in a manner which depends upon 
the magnetization characteristic. The range of forces over which our 
equations are valid depends simply upon the number of terms taken 
in the development of the branch equations, and is not necessarily 
restricted to small forces. The hysteresis loop coefficient enters into 
the impedance offered to the fundamental frequency by a coil enclosing 
the core material in question in such a way that the third harmonic 
produced may be deduced from the change of resistance with current. 
Further, the ratio of the change of reactance with current to the 
change of resistance with current may be used to provide a test of 
Rayleigh's relation, which is found to hold for some materials, while 
it is invalid for others. The precision obtainable in the evaluation of 
hysteresis loop coefficients is much greater by the a.c. bridge measure- 
ment under proper experimental conditions than by the analysis of 
ballistic loops. Incidentally attempts have been made to obtain B-H 
loops by AC methods with the aid of the Braun tube, for example, 11 
but the precision attainable is not sufficiently high for our purpose. 

Complex Magnetizing Force. Our preceding analysis has furnished 
us with the fundamental voltage drop across a coil enclosing the iron 
core under consideration, together with the third harmonic voltage 
generated in the coil winding due in general partly to the non-linear 
B-II relation and partly to the effect of hysteresis. The generated 
voltages corresponding to the fifth, seventh, and higher orders are 
also calculable by the same methods but will not be specifically con- 
sidered since they are smaller than the third harmonic at low fields. 
This generated third harmonic voltage exists in its entirety across the 
coil winding only when the impedance of the external circuit is much 
higher than that of the coil at the harmonic frequency, a condition 
not usually satisfied in telephone circuits. A current of the third 
harmonic frequency then flows in the circuit, its amplitude and phase 
depending evidently upon the generated third harmonic voltage — 
that is, the coil structure and core material — as well as upon the total 

11 Peterson, Phys. Rev., Vol. 27, No. 3, p. 320. 
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circuit impedance to the third harmonic which includes that of the coil. 
We have now to determine the coil impedance to a complex magnet- 
izing force constituted by the fundamental and the third harmonic. 
It is evident at the outset that if we restrict consideration to a very 
small third harmonic, the impedance to the fundamental cannot be 
very materially altered. 

In general the fundamental and third harmonic may exist in any 
phase but for our present purpose we shall take the magnetizing force 
to be 

H = Hi cos pt + H 3 cos 3pt, (29) 

in which the phase angle is assumed zero. This seems to be an 
arbitrary assumption, but it may be shown that the results are not 
materially different for other phase displacements, and there is some 
gain toward simplicity of treatment by taking the phase angle zero. 
The equations of Part 1 may then be carried over without change, and 
details of the analysis are given in Appendix 4. 

The ratio of the resistances to the third harmonic and to the funda- 
mental is 0.77, from the relations deduced in the appendix. Some 
experimental work has been done to test the validity of the expres- 
sions derived above for the reactance and resistance components of 
the third harmonic, the results of which are given in Fig. 4. It is seen 
that the check for the inductance to the third harmonic is very good 
but that the resistance component appears to be substantially that of 
the fundamental rather than 77 per cent of it, as the analysis indicates. 
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Fig. 45 

The experimental method for obtaining this result — the impedance to 
a small third harmonic in the presence of a relatively large funda- 
mental — is illustrated in Fig. 5. The method consists in measuring 
the third harmonic current by means of a current analyzer 12 for a 
number of circuit conditions in which the fundamental current is 



TO CURRENT 
ANALYZER 



X HIGH £ OUTSIDE BAND 
Fig. 5 

maintained constant. The circuit is first tuned to the third harmonic 
by varying the capacity C in the third harmonic circuit, and the current 
is then measured for a series of values of the series resistance r. A 
shunt resonant circuit tuned to the fundamental is inserted in the 
third harmonic path so as to separate effectively the third harmonic 
circuit from that of the fundamental. With this precaution taken to 
avoid harmonic production in the analyzer and to maintain the funda- 
mental current constant while r and C are varied, the inductance 
to the third harmonic is obtained from the resonating capacity, and the 
resistance is determined as that value of r for which the third harmonic 
current falls to half its maximum. 

12 For details of current analysis see the paper by A. G. Landeen, B. S. T. J., 
April 1927. 
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Part 3. Application of the Analysis 

Air-Gaps and Dilution. Under certain conditions improvement in 
the operation of iron core coils toward freedom from harmonic pro- 
duction may be attained by inserting air-gaps in the magnetic path. 
The expressions which we have derived up to this point are valid for 
a material having the constants assigned, and the question now arises 
as to the parameters which characterize the operation of the iron core 
including air-gaps, and their relation to the parameters for the original 
core without air-gaps. With these relations given, our previous work 
may be applied to cores with air-gaps. 

In establishing the correspondence between the parameters for the 
two cases, it is instructive to use two methods — one a direct attack, 1 '' 
the other resting on an analogy with non-linear vacuum tube circuits. 14 
We may determine the effects sought for by the direct method on 
consideration of a single branch of a hysteresis loop, which is expressed 
by equations (4a) or (5a) of Part 1, 

B(h, II) = EJjhJrB- ( 36 ) 

Now with an air-gap in the magnetic circuit, the magnetomotive force 
effective is that applied, reduced by the drop across the air-gap, or 

m' = m — p<p, /,«■> 

M' = M - p$, K J 

in which mM, <p$ are instantaneous and maximum values, respec- 
tively, of the impressed m.m.f . and flux, and in which p represents the 
air-gap reluctance 

P = \/A, (38) 

X being the length of air-gap and A the core cross-section. 

In order to apply (37) we re-express (36) in terms of magneto- 
motive force and flux as follows : 



A a, 



<p(m,M) = EE^m'T, 

r s * 



(39) 



where I is the length of magnetic circuit in the iron. Equations (37) 
may now be substituted in (39) to yield 



,Aa 



<p(m, M) = EZ T=? (™ ~ WW*" - P*Y- ( 4 °) 



T S 



/'+• 



13 Due to Mr. H. P. Evans. 

14 See Appendix 5. 
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This equation may be solved for tp by identifying coefficients when we 
substitute the solution 

rill. M) = ZZ^f™ r M>, 

T S * 

which is equivalent to 

B(h, H) = ZJjh.'hrtP. (41) 

r s 

Carrying out the operations, the primed coefficients are determined 
in terms of the original coefficients and the core structure as follows: 

flio' = flio /( 1 +7 a io ). 

a 2 = — a 20 ' = co2 / f 1 + j- aio J , (42) 

a n ' = On/ ( 1 +|-OioJ • 

It is clear that the effect of an air-gap is to diminish the higher 
order coefficients to a greater extent than those of lower order; no 
new coefficients are introduced. If we refer to a constant flux density, 
then the impressed force is 

£r-J5r(l+£a 10 ) (43) 

and the modulation voltage becomes proportional to 

a 02 'H' 2 =- tv// " J 



1 _!_ X 
1 +T a io 



so that the harmonic e.m.f. has been reduced in the ratio of 1 + j am 

to unity, which is precisely the ratio of initial permeabilities. 

Dilution. If the magnetic material is diluted with a non-magnetic 
substance (insulation) so that the effective air-gap length is X, the 
above equations for the flux density still hold. The effective cross- 
sectional area, however, is reduced to 

A' = A^i^J, (44) 

which is ordinarily of small account compared to the other factors 
involved. 

Choice of Core Material. The ideal characteristics for a core ma- 
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terial in respect to its freedom from harmonic production may be 
determined by a consideration of the circuit problem. The third 
harmonic driving e.m.f. is obtained from the last section as 

E 3 = pnAa 02 H*10-*, 

^n*A (45) 

in which n represents the number of turns enclosing the core of area 
A, d is the mean diameter of the toroidal core and a 02 is the hysteresis 
coefficient. 

Suppose that we have to design a coil of fixed inductance in which 
the harmonic is to be a minimum, so that 

L = K 7 ^ = const. (46) 

We proceed to the consideration of a number of special cases subject 
to condition (46). 

Case 1 — L Fixed, n Variable. From (46) 

n = const. — r I , 

which may be substituted in (45) to give 

Hence in a coil of fixed inductance, fixed core area, and fixed core 
diameter, minimum harmonic voltage is produced with a material 
for which a 2/V' 2 is minimum. 

Case 2 — L Fixed, A Variable. From (46) 

A = const, d/n 2 n, 

which, upon insertion in (45), yields 

in which the modulated voltage varies linearly with a 2//J- 

Case 3 — L Fixed, d Variable. By a similar procedure, we obtain an 
expression for the generated third harmonic voltage 

For further work we suppose the coil reactance considerably greater 
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than the circuit resistance at the fundamental frequency, and consider 
the case of an input transformer in which the secondary is practically 
open-circuited. 

Case 4 — Input Transformer. In accordance with the above assump- 
tion the fundamental current through the coil is determined by the 
coil reactance x or 

d 



T - £ - 

* 1 ■> A ' 



Putting this in (45), we find 



r, . 1 fl 02 

1lA IT 



(50) 



which is identical with Equation (49). 

In the four cases treated above it has appeared that there are three 
quantities which characterize the modulating properties of a ferromag- 
netic coil under different conditions — a^ln, a 02 /tJ. 312 , aoi/fj?. The values 
of these ratios have been determined for a number of materials and are 
tabulated below. 



002/M 



(102/n 



3/2 



a ul /n" 



B dust 

C dust 

Permalloy dust ,5 
Perminvar 15 . . . . 
Silicon Steel .... 



18 X 10- 3 
11 
17 
2.1 
2500 



3.1 X 10- 
2.2 

2.1 
0.1 
210 



0.52 X 10- 
0.42 
0.24 
0.0045 
16 



These results apply only when eddy currents are negligible; as the 
frequency is raised eddy currents effectively screen the core and the 
harmonic flux is reduced to a greater extent than is the fundamental. 
This effect, as is well known, depends upon the specific resistivity of 
the core material and will not be evaluated here. 

It has been pointed out that under different circuit conditions an 
index of the modulation is given by the ratios a^Jn, ao2,V /2 . aoi/tr, 
depending upon specific conditions. For diluted materials, these ratios 
referred to the undiluted material become 



002 



a 2 



Q02 



010(1+7-010) a,o(l+7;aioj flio 2 ( 1 +^-aio j 

and the utility of air-gaps toward the reduction of harmonic is made 
evident quantitatively in every case. 

In order to test these relations a comparison was made of the 

" Laboratory specimens. 
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coefficients for specimens of grade "B" and grade "C" dust which 
represent different dilutions of electrolytic iron. The coefficients 
used were those tabulated in Part 1. In accordance with equations 
(42) above we may write 

flio/(l+jflio) =26.2, a 10 /(l +ja l0 ) =35.3, 

an/(l+ja lo y= 0.59, a„ /( 1 + j a 1Q J = 1.53, (51) 

a m IU +jflio) 3 = 0.29, fl «/( 1 +T ai0 ) 3 = °- 65 - 

From these, we should have equality of the ratios 

0.59/1.53, 0.29/0.65, (26.2/35.3) 3 , 
or 

0.37, 0.45, 0.41, 

which are evidently in fair agreement. It is clear then that the 
properties of diluted materials may be calculated from the character- 
istics of the original material, at least when the process of dilution 
does not change the intrinsic properties of the magnetic material 
involved. 

Applications to Transformers. The third harmonic flux component 
may be obtained when the fundamental magnetizing force is given, 
and the latter is easily obtained in toroidal core inductance coils from 
the relation 

h = OAni/d, (52) 

where both h and i refer to the fundamental frequency. In a trans- 
former, however, the net magnetizing force is obtained as the sum of 
two components, one due to the primary and the other produced by 
the secondary, so that some further investigation is required before 
the net magnetizing force in the core may be calculated in terms of the 
transformer constants and the primary current. 

Net Magnetizing Force. To obtain the net magnetizing force which 
determines uniquely the generated flux components we are required to 
obtain the primary and secondary currents (ifc), to multiply each 
current by the number of times it encircles the core, to add the two 
products, and finally to multiply the sum by 0A/d: 

flnet = 0.4(«i«i + n 2 i 2 )/d. (53) 

To evaluate (53) then we shall solve for the primary and secondary 
fundamental currents in terms of the circuit constants and the applied 
potential. 
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The transformer circuit equations from which the currents may be 
evaluated are the familiar ones 

E = Ziii + jpMh, ~ 4 * 

= Z 2 i 2 + jpMi u 

in which we assume the transformer to be an idealized structure — the 
capacity and leakage effects of actual transformers will be assumed 
combined with the connected impedances for simplicity — they will 
therefore not be dealt with explicitly. From the second of (54) is 
obtained the relation between the two currents 

*a = - jpMii/Zt (55) 

which may be substituted in the first of (54) to give the usual expression 
for the primary current 



;. - » / 



i\ = E 



(56) 



I Ri + -^- r K 2 -tJ I A] YT 

An expression for the net magnetizing force in terms of ii may be 
had putting (55) and (56) in (53) 

HM = oMA^_jp^y (57) 

in which K represents the turns ratio : 

K=^ = ^l- (58) 

«i -4x1 



If we make the convenient assumption, which is closely approximated 
under working conditions, that the coupling is perfect (no leakage) we 

have 

pM = -\XxXo (59) 

and (57) may be simplified to the form 

^=9^*2. (60) 

a Z.2 

This equation shows that the net magnetizing force may be obtained 
from that of the primary alone by applying the reduction factor R 2 /Z 2 . 
Now in well designed transformers the winding reactances are much 
greater than the resistances to which they are connected, and in this 

case 

X\ 2 »R\ 2 . (61) 

Applying this further simplification to (60) we have finally for the net 
50 
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magnetizing force 

OAtiiii R 2 . . 

-«net = — 2 — X~ ' * ' 

in which the reduction factor applied to the force due to the primary 
alone is seen to be the ratio of resistance to reactance in the secondary 
circuit. 

This equation with the aid of (56) may now be put in terms of the 
primary voltage — a fixed quantity — rather than in terms of the 
primary current which is variable according to the circuit resistances 
used. Applying the assumptions (59) and (61), Equation (56) may 
be written 

n = EliR, + R 2 /K 2 ), (56a 

which may be put in (60a) to express the net force explicitly in terms 
of the circuit constants: 

_ 0.4w,E / 1 \ 

i/ "" t dX l \ 1 + KtRy/Rt J {0Z) 

In view of (58). When the transformer is terminated in its normal 
resistances the expression within parentheses reduces to the value one- 
half. It is of interest in this connection to compare the field here with 
that produced with the secondary open ; in the latter case it is clearly 

OAniE/dXi, 

which is twice as great as the field in the properly terminated trans- 
former. This result is otherwise evident, for in the properly termi- 
nated transformer but half the generator voltage is applied across 
the primary. 

Third Harmonic. According to Equation (28) the amplitude of 
third harmonic voltage generated in a coil of n turns encircling a core 
subjected to a magnetizing force of H net gilberts/cm is 

E 3 =^- lO-*a 02 pnAH 2 net (63) 

JTT 

and the generated primary or secondary voltage is obtained by attach- 
ing appropriate subscripts to n. Because of the coupling, the two 
circuits react on one another at the third harmonic frequency as they 
did at the fundamental frequency, and to find the two third harmonic 
currents we must go through a procedure analogous to that followed 
for the fundamentals except that here we have an independent 
generator voltage effective in each of the two windings. 

Indicating quantities referring to the third harmonic frequency 
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by the subscripts p, s for primary and secondary circuits, respectively, 
we have the circuit equations for the third harmonic currents: 

e p = Z p i p + j3pMi s , 
Ke p = Z.i, +j3pMi p , K ' 

and by equating the two expressions for e p we obtain a relation between 

the two currents : 

. KZ p -j3pM. ,,, x 

*' = Z.-jipMK ^ (65 > 

Putting i p in terms of e p , and using (59) and (61) we find 

i, = e p KR P IR t X p {\ + K 2 R pl 'R t ). (66) 

Suppose now that we are dealing with pure resistance loads so that 

Ri = R P , i?2 = Rf 

We may then substitute (63) for e p in (66) to obtain the third harmonic 
current in the secondary. 
From (63) 

e P - ^ 10 "anPmA -fzj- { - t + R2Ri/Ri j • (67) 

and by substitution in (66) 

E*KR ia 

•• = 7 &\* ' < 68 ) 

SXJUl 1 +JP^ 

If we consider the factor containing the resistances, 

R 2 (l + K 2 i?!/i? 2 ) 3 ^° v; 

it is zero for Ri/Ro zero or infinite, and reaches a maximum under the 
condition 

mm - °- ™ 

The third harmonic secondary current is maximum when 

Now i£ 2 is fixed by the turns ratio (58), so that we may say the 
secondary third harmonic current is maximum when the primary 
resistance is made half its nominal value, or when the secondary 
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resistance is made twice its nominal value. This is well borne out by 
Fig. 6 due to A. G. Landeen taken under conditions to which the 
above discussion applies. 
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Fig. 6 

At first sight the result obtained, in which the secondary harmonic 
current increases as the secondary resistance is increased for all values 
of secondary resistance below twice the nominal secondary resistance, 
seems a bit strange since the fundamental secondary current decreased 
under the same conditions. A little thought however shows that the 
increase of secondary resistance acts in two ways ; first to reduce the 
harmonic current directly, and second to increase the net magnetizing 
force in the core which in turn increases the generated harmonic. 
Thus with zero secondary resistance the net magnetizing force is zero 
and no harmonic flux can be produced, while with large secondary 
loads the flux density in the core or the net magnetizing force reaches 
an asymptotic maximum so that the generated third harmonic increases 
slowly if at all while the reduction in secondary current by the re- 
sistance is continuously increasing. Under the assumptions noted 
Equation (68) is valid for any circuit condition and it will be observed 
that i, decreases with the primary resistance below the optimum point. 
The value of a is explicitly 

ao2 ( fii 3 A 



a = 2.72 X 10- 10 



p 2 LS 



d 2 



(72) 
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in which the bracketed expression coincides with the variable factor 
in equation (45). Hence when L\ is constant the choice of core 
material follows the rules laid down for inductance coils. When 
the inductance is not constant, however, the factor becomes pro- 
portional to 

a = const. 3 °' (73) 

with the understanding that the turns ratio, resistances and frequency 
are fixed. Thus the secondary harmonic current with a given material 
is reduced by increasing the turns and core area and reducing the 
diameter. As far as core materials go, we have for the significant ratio 
«02/m 3 the values: 

Material a i/n 3 

Silicon Steel 130 X 10"* 

Permalloy Dust 16 2.4 X 10"* 

B Dust 5.3 X 10"* 

C Dust 4.3 X 10" 1 

Perminvar » 0.05 X 10^ 

The great superiority of perminvar indicated above is restricted of 
course to fields of the order of less than 0.1 gilbert/cm. above which 
it tends to become smaller; at a field of 0.7 for example the above 
factor for perminvar would be multiplied by the factor three. Perm- 
alloy dust is observed to be approximately twice as good as the iron 
dust cores. 

A very important question to answer with regard to transformer 
cores is this — what benefit can be gained as to harmonic production by 
inserting air-gaps, or by diluting the core material, leaving everything 
else unchanged. This is evidently to be answered by an investigation 
of the ratios a 2/M 3 and a' 02 /n' 3 , the primes referring as before to the 
diluted material. These relations are given by Equation (42) in 
which n and a i0 are interchangeable: 



aio' = Cio / f 1 + y aio j = /*'» 
002 = flo2 / I 1 + 7 aio I • 



These permit us to evaluate a'oo/ju' •" 



1 _L X V 

ao2 a 2 V ' / a ° 2 (7i\ 

M / 1 J- X n \ ai ° M 

I 1 + j aio 



16 Laboratory specimens. 
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which clearly indicates that no change is produced in the secondary 
third harmonic current by introducing air-gaps or by diluting the core 
material, the change in core area being negligible. This relation has 
been verified experimentally. 

The discussion above considers a constant potential applied to a 
transformer circuit and the relations derived are somewhat changed 
when we consider constant current or constant potential transformers. 
We seldom have to do with constant current transformers but constant 
potential transformers are met with frequently in vacuum tube circuits 
as input transformers or interstage transformers. The above discus- 
sion may be applied to this case by taking the primary generator 
resistance much lower than its nominal value — a condition, inci- 
dentally, which works toward the suppression of harmonic. 

The conclusions are also somewhat altered when we have networks 
offering different impedances to the harmonic and the fundamental. 
Thus it is evident that the third harmonic current can be eliminated 
from both primary and secondary circuits by inserting a high series im- 
pedance to the harmonic frequency, or by encircling the core with a 
winding connected to a network which has a very low impedance to the 
third harmonic and high impedance to the fundamental. Further, in 
single frequency transmission the result may be equally well obtained 
by shunting a series tuned circuit around the primary or secondary 
winding. 

Applications to Harmonic Production in Loaded Lines 

If distortion takes place at any one point of a loaded line, the amount 
of distorted current received at the far end depends upon the ampli- 
tudes of the currents producing it, and upon the line attenuation from 
the point of origin to the far end. The phase similarly depends upon 
the phase of the fundamentals, and upon the phase shift of the line. 
When we have a number of distorting sources at different points along 
the line, the net distorted output is obtained by combining vectorially 
the currents due to the individual centers of distortion, since it may 
be assumed that no interaction exists. In a uniformly loaded line 
we may think of these sources of distortion as being uniformly dis- 
tributed, but the amount of distortion introduced is, on the contrary, 
not distributed uniformly on account of the line attenuation which 
reduces the distortion generated at the far end of the line. 

It is apparent, then, that if we are to calculate the net distortion 
introduced by the line, a complete specification of the phase shift 
and attenuation is required, together with a knowledge of the law of 
production of the distortion. This last also requires a specification of 
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amplitude and phase angle in their relation to the fundamentals 
producing the distortion. 

A situation somewhat analogous to the one under discussion which 
is less involved, however, is found in the diffraction of light, in which 
the illumination at a point proceeds from a number of coherent sources 
at various distances from the point. An important difference in the 
two cases is the attenuation of the transmitting medium which is 
ordinarily small in the optical case, so that the results of the optical 
investigations cannot be taken over directly. 

In the following we propose to calculate the third harmonic output 
currents with the aid of our previously established relations. The law 
of harmonic production has been quite definitely established in the 
low frequency-low flux density range. The driving third harmonic 
voltage of frequency 3/ produced by a fundamental current of rms 
value 7 is given by (25) and (26a) as 

«M - — 

E 3 = 0.72 X 10- 8 a 02 ^// 2 - MI *> ( 75 ) 

while the phase angle of the third harmonic generated potential is 
related to that of the fundamental current by the expression 

03 = 30i. (76) 

The above data are sufficient for a solution of the problem of distortion 
in continuously loaded lines when the propagation constant is known 
as a function of frequency. 

The fundamental current at any point distant x from the sending 
end of the continuously loaded, properly terminated line is 

/ = Joe-'*, (77) 

where 

P = a +7/3. 

The third harmonic driving e.m.f. dE 3 generated in a length dx of 
the line may be written with the aid of (75), (76), and (77) as 

dE 3 = M/ 2 e- (2Qr+ ^ )x . (78) 

Writing the line impedance as Z , the resulting current at x is 

dl 3 =Mll e -v«+m*. (79) 

With the line parameters a function of frequency we may write for the 
propagation of third harmonic current 

i z = I 3 t-(y+i»*. (80) 
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The distortion at x produces a third harmonic current at the receiving 
end 

2Z 



di t = I ^±L e -^+flfiu € -iy+im-x) t ( 81 ) 



The total third harmonic current at the receiving end may now be 
obtained by integrating (81) over the length of the line /, which gives us 

2Z [(2a - 7) +j(3I3 - 8) J 

Inasmuch as we are concerned with the output amplitude, the above 
expression may be put in a somewhat more convenient form: 

I- 12 = (MUX. e ~ 2yl 

l * 31 ^ 2Z ) (2a - 7 ) 2 + (3/3 - 5) 2 

X [1 + e -'^ 2a - y)l — 2e-< 2 ?-T» cos (3/3 - 6)1]. (83) 

Thus when I is zero the harmonic vanishes as it should, and as / 
increases the current passes through maxima and minima determined 
by the cosine term. If the attenuation is not very great, the maxima 
occur approximately at the line lengths 

_ (2n - 1)tt 
3/3 - 5 ' 

where n is a positive integer. These distances correspond to odd half 
wave-lengths, as is true of the optical case. As I is increased the 
bracketed expression approaches unity and the current falls expo- 
nentially. Before this point is reached, however, the current increases 
in certain regions as the line length is increased. The fact that in an 
actual line the parameters vary with frequency means that these 
maxima and minima will vary in position according to the frequency, 
so that a maximum for one frequency may well coincide with a mini- 
mum for another frequency. 

In the case of lumped loading, the integrations used for continuous 
loading are replaced by summations, as was done by Mason in the 
unpublished investigation previously cited. 4 If the spacing between 
coils is Xi, and if we have n coils for which nx\ = y, the received third 
harmonic current at the end of a properly terminated line may be 
written 

. _ MIp 2 e-^+' i)l {l - c-Efa-T+yPg-mv} 
H ~ 2Z (^-r+jW-Wx, _ 1) ' t°4) 

which is to be compared with (82) for the continuously loaded line. 
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Appendix 1: Simplification of Loop Equations 

From the first of the three properties mentioned we have a 00 = 0, 
and from the second property 

B x (h,H) = - B 2 (-h,H) (3) 

whence, if we write 

Bi(k, H) = a 10 h 4- a 01 H 

+ a 2Q h 2 + a n hll + a 02 -f/ 2 

+ a 30 h 3 + a 2i h 2 H + anhH 2 + a 03 H 3 • • • , (4) 
then 

B 2 (h, H) = aioh — aoiH 

— a 20 h 2 + anhH — a 02 H 2 

+ a 30 h 3 - OiiPH + anhli 2 - a 03 H 3 • • • . (5) 

From the third property, the two branches meet at the loop tip which 
lies on the magnetization curve for which h = H, or 

3i(ff, H) = B 2 (H, H). (6) 

From the two equations (4) and (5) we have by virtue of this relation 

a 01 H + (a 20 + a 02 )H 2 + (a 21 + a 03 )H 3 = 

and since this relation holds for every value of II, the coefficients of 
each power of H must be zero so that 

floi = = aoo, 

002 = — #20, (' ) 

fl03 = — #21- 

The final expressions for the branches are then evidently obtained by 
putting (7) in (4) and (5) which become 

Bi(h, H) = aioh 

— a 02 h 2 + a u hH + a 02 IP 

+ a 30 h 3 - a 03 h 2 H + a l2 hH 2 + a 03 # 3 , (4a) 

J5 2 (/t, 77) = auh 

+ a 2^ 2 + anhH — a 02 H 2 

+ a 30 /* 3 + a 03 h 2 II + ai 2 MF - a 03 H 3 . (5a) 

Appendix 2: Rayleigh's Relation 

If we suppose the lower branch of any loop, when referred to the tip 
of the largest loop considered, to be given by the equation 

B x = mo/*! + vhi 2 + \h? + ah? (12) 
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we may refer the family of branches to the origin by the transformation 

I : t 1 - f ■ d3) 

if B m II refer the midpoint of the largest loop to the tip. Then 

B m = i* H + 2vIP + 4X# 3 + 8a,H 4 . (14) 

Putting (13) in (12) 

B +B n = n (h + II) + v(h + tf) 2 + X(A + i?) 3 + «(A + II) 4 , 

whence, subtracting (14), 

B' = »oh 4- v(h? + 2hH - H 2 ) + \(h 3 + 3h 2 H + 3/*H 2 - 3H 3 ) 

+ co(/? 4 + 4MH" + 6WP + 4/z# 3 - 7 II 4 ), (15) 

which represents the hysteresis branch equation referred to the origin, 
on the basis of loop similarity. 

The coefficients obtained by the two methods may now be compared. 
Thus identifying coefficients of (15) with those of the general equation 

(1) 

010 = Mo. #n = 2v, #02 = — U O12 = 3X, 

a 2 o = v, a 2i = 3X, a 3 = 3X, a i3 = 4co, (16) 

#30 = X, #31 = 4t0, fl04 = "- 7<o, #22 = 6a>. 

a4o = co, 

Appendix 3: Alternating Magnetization, Sinusoidal 
Magnetizing Force 

The resulting expression is simplified if we make the following 
substitutions 

a = a 02 H 2 + a 03 H 3 = B(0, II), 

= oio + a n H + anil 2 , (19) 

5 = a 30 tf 3 . 

It may be noted that a is the remanence and that /3 is an approxima- 
tion to the permeability, in fact the permeability is given as the sum 
of /3 and 5. With (18) and (19) inserted in the branch equations, then, 
we have 

Bi(H cos pt, H) = a + /3 cos pt — a cos 2 pt + 5 cos 3 pt, 
B 2 (H cos pt, II) = - a + /3 cos pt + a cos 2 £/ + 5 cos 3 ££. 

For convenience we shall express these relations in terms of multiple 
angles, and we have for the equation of the upper loop family 

Bi(Hcos pt, H) = - a/2 + (/3 + 35/4) cos pt - a/2 cos 2pt + 5/4 cos 3pt. 
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If we write 

= a/2 = (fl.02// 2 + a 03 H 3 )/2, 

(20) 



A = a/2 = (a 02 // 2 + a 03 H*)/2, 

B = + 35/4 = a 10 + auff + a 12 fl* + 3a ao H*/4, 

C = = - A, 

D = 5/4 = a 3 o#74, 



the final form for the loop equations is 

Bi(H cos pi, II) = A + B 17 cos /»/ + C cos 2# + D cos 3pt, , . 
5 2 (#" cos pi, II) = - A + B cos pt - C cos 2/>* + D cos 3#. ^ ' 

We are now in position to combine the two equations of (21) in a 
Fourier series valid over the entire cycle as 



B = j + £(&* cos kpl + a k sin kpt), (22) 



where 



1 f 2r 
0* = - I /(/>0 sin kpt d(pt), 

■" Jo 

&* = ~ f 2 '/(^) 



(22a) 



cos £/>/ dipt). 



For our particular case we have, since B = Bi for the first half of the 
cycle, and B = B 2 for the second : 

Trafc = I 5 a (fc, #) sin kpt d{pt) + I Bi(h, H) sin kpt dipt), 

J-T JO 

r&i = I 5 2 (//, II) cos jfe£/ <*(£/) + I 5i(A, i7) cos ^ dipt). 

These integrals may be simplified considerably when we take advantage 
of the fact that both Bi and B 2 are even functions of the time as given 
by (21). Thus 

■Ka k = j [Bi(h, II) - B 2 (h, H)2 sin kpt dipt), 

Jo 

irb k = f lBi(h, H) + B t (k, H)] cos kpt dipt). 

Jo 
Referring to (21) we may then write 

a k = - I (A + C cos Ipt) sin kpt dipt), 

rJo (23) 

2 C 

b h = - (5 cos pt + D cos 3pt) cos kpt dipt). 

TT Jo 

17 This coefficient is not to be confounded with the general expression for flux 
density. 
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Upon integration of (23) the coefficients for the fundamental and third 
harmonic flux components are found to be as follows : 

ai = -(A - C/3), 6i = 5, 

(24) 

which, by reference to (20), may be put in terms of the branch coeffi- 
cients. 

Appendix 4 — Impedance Reaction to a Small Third Harmonic 
in the Presence of a Large Fundamental 

We have for the two hysteresis branch equations from Eqs. (4a), 
(5a) 

Bi(h, H) = B(0, II) + /3/» + yh* + a 30 h 3 + • • •, 
B 2 {h, H) = - B(0, H) + fih - yh* + a 30 h 3 + • • •, 



(30) 



(31) 



in which 

= aio + auH -f- a n H 2 , 
7 = — (ao2 + a<tzH)- 

Putting (29) in (30) we get 

B(h, H) = A + B cos pi + C cos 2pt + 2) cos 3# + F cos «/>* 
+ G[cos (n + l)/>/ + cos (n - l)pt] 

+ /[cos (n + 2)# + cos (» - 2)#] (32) 

in which the coefficients have the following significance 

A = B{0, II) + .fti 2 7/2, F - pH 3 + 3a i0 Hi 2 H a i2, 

5 = jStfi + 3a 30 tfi 3 /4, G = 7^1-^3, / 33) 

C = 7#i 2 /2, 7 = 3a 30 tfi 2 // 3 /4. 

D = a 30 -H"i 3 /4, 

The coefficients of the Fourier Series for the output wave may 
now be obtained as before by combining the two equations (30) since 
each one is operative during one-half the cycle. There results an 
expression similar to the one obtained in the single frequency case, 
and since we have 

/ = b /2 + 2a* sin kpt + 26* cos kpt 
the coefficients are evaluated from the expressions 
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? r 2r 

a k = -\ {A + C cos Ipl + G(cos Apt + cos 2pt)) sin kpt dipt), 

■JTJO 

b h = - I (-B cos p/ + (D + F) cos 3£/ 

ttJo 

+ /(cos 5^>/ + cos />/)) cos kpt d(pt). 

Upon integration we find 

6, = 5, & 3 = F. (35) 

Comparing these two coefficients we see that at low amplitudes we 
may write 

bi = iiH\, b 3 = nH 3 , 

in which the permeability is the same to the two components, and is 
determined by the fundamental amplitude. For the dissipative 
terms we find 

ai = -{A- C/3 - 2G/5), 

* (35') 

a» = - (A/3 - 3C/5 + 6G/35), 

7T 

but some care is required in interpreting these expressions. Inasmuch 
as we are primarily interested here in determining the dissipative 
component to a third harmonic magnetizing force of amplitude TI 3 , 
we are required to select from a 3 only those terms containing H 3 , 
which means the single term 24G/357r. The other terms take care 
of the harmonic producing properties of the core and do not affect 
the impedance to the third harmonic. The impedance term for the 
third harmonic comes down to 

-^-(a 2niH 3 + a 03 IISH 3 ), 

307T 



which may be written as 



24 B(0, H) 
35r ' Hi 



This may be compared with the corresponding term for the funda- 
mental given by (25). 

Appendix 5. Effect of Air-Gap by Vacuum Tube Analogy 

In the elementary treatment of non-linear two element vacuum 
tube circuits, approximate solutions are obtained in the form 

/ = /l + Jt + • • • J; 
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where / represents the variable part of the space current on the basis 
that J n represents the rath approximation, but that the series converges 
rapidly so that we need consider only the first two terms to arrive at 
a substantially accurate result. The expressions derived for the first 
and second approximations are known to be 

/, = a x EliX + ai-R). 
J 2 = a 2 E*/(l + ai R)», 

where the 'a' coefficients describe the tube characteristic 

J = ay -f- a 2 v 2 , 

R being the external plate circuit resistance, v the tube potential, 
and E the circuit e.m.f. 

Turning now to the equations for the hysteresis loop branches, we 
have from (4a) 

B = a w h — a 02 k 2 + a u hH + a^H 2 . 

Hence by the analogy between flux and current, and between reluc- 
tance and resistance, the first order terms are reduced by the factor 
1/(1 + a-iR) which corresponds to 1/(1 + \ai /A), and the second 
order terms are reduced by the cube of this factor, which yields the 
same results (42) as the laborious direct method. 



